Abstract 


This  report  describes  the  development  of  an  alternate  formulation  of  the  QD 
digital  filter  using  matrix-vector  notation  and  methods.  The  objectives  of  the 
report  are  to  clarify  the  concepts  on  which'  the  filter  is  based  and  to  briefly 
indicate  the  structural  similarity  of  the  QD  and  Kalman  filters  as  well  as  the 
difference  between  them.  The  technique  for  inserting  into  QD  observed  derivative 
data  obtained  independently  of  the  normal  input  observations  is  briefly  described. 
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A  Matrix-Vector  Formulation  of  the  QD  Digital  Filter 


Introduction 

In  July  1967  the  writer  developed  a  new,  more  general  formulation  of  a 
very  efficient  type  of  real-time  digital  filter  and  called'  it  the  "QD  Filter". 
Actually,  the  structures  of  the  computational  formulas  for  the  first  and  second 
order  QD  filters  were  not  new.  As  early  as  1957,(2)0)  for  example,  workers  in 
radar  engineering  developed  the  well-known  a-$  digital  filter  which  is  used 
primarily  to  predict  estimates  of  position  coordinates.  The  ot-6  filter  is  com¬ 
parable  to  a  first-order  QD  type  of  filter.  Then  in  1964  J.  J.  Lynn  O)  developed 
a  second-order  filter  with  three  coefficients  a,$,  and  Y  , analytically  determined 
to  attain  the  optimum  trade-off  between  distortion  of  the  true  input  information 
and  attenuation  of  the  corrupting  noise  inherent  in  the  input  data.  The  QD 
formulation,  in  contrast  to  these  earlier  developments,  is  jiuch  more  general 
because  it  is  based  on  ap  N-th  order  polynominal  fit  of  a  span  of  input  data  points, 
subject  to  the  least-squares  error  criterion  and  intercept  and  slope  constraints. 
Exhaustive  experiments  have  shown  that  the  performances  of  the  second  and  third 
order  QD  filters  are  essentially  the  same  as  the  respective5 performances  of  the 
equivalent  classical  constrained  least-squares  filters. 

The  minimal  computing  time  required  by  the  QD  filter  is  independent  of  the 
point  span  length.  Because  of  this  computational  efficiency,  the  second-order 
QD  filter  is  now  implemented  in  a  variety  of  computer  programs  for  filtering 
radar  data  at  White  Sands  Missile  Range,  e.g.,  the  real-time  program  for  flight- 
safety  support  of  PERSHING  missile  overflights.  The  third-order  version  has 
been  used  experimentally  and  displays  negligible' "velocity  lag"  error  in  velocity 
estimates  obtained  from  radars  observing  a  missile. during  its  thrusting  periods. 

It  is  understood  that  Lynn’s  cx,$,Y  digital  filter  has  been  used  in  similiar  real¬ 
time  computer  programs  at  the  Eastern  Test  Range. 

It  has  been  observed  that  the  QD  computational  formulas  have  the  same  general 
structure  as  those  of  the  discrete  Kalman  Filter.  This  similarity  is  quite  natural 
because  both  filters  are  recursive.  Otherwise,  they  are  so  different  that  the  QD 
filter  can  not  be  considered  as  a  special  case  of  the  Kalman  filter,  which  will  be 
quite  apparent  to  readers  who  are  familiar  with  the  Kalman  theory. 

The  original  QD  filter  development  employed  classical  summation  (scalar) 
notation  which  was  extremely  tedious  for  the  reader  to  follow.  It  is  hoped  that 
this  new  development  using  matrix-vector  notation  presented  in  the  following 
chapters  will  be  clearer  and  facilitate  much  easier  reading.  The  manner  in  which 
the  development  is  presented  has  several  other  specific  objectives:  (1)  to  show 
the  Kaiman-knowledgable  readers  that  the  QD  filter  can  be  cast  in  matrix-vector 
notation  and  that,  despite  the  structural  similarity,  the  QD  theory  is  distinctly 
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Notation 


In  our  development  we  will  use  the  following  notation: 

1  -  Matrices  In  capital  letters  cr  matrix  algebraic  expressions  in 

brackets; 

2  -  Column  vectors  in  lower-case  letters; 

3  -  Row  vectors  in  lower  case  letters  with  a  "T"  superscript; 

4  -  Scalars,  vector  elements,  matrix  elements  in  lower  case  Greek  letters; 

5  -  Minus  one  superscript  for  the  inverse  of  a  square  matrix; 

6  -  Subscripts  on  derivative  vectors  denote  the  point  in  the  span  at  which 

the  derivatives  are  evaluated; 

7  -  Subscripts  will  also  be  used  to  designate  the  components  of  partitioned 

matri'ces  and  vectors; 

8  -  Estimates  will  be  designated  with  a  "hat"  and  predicted  values  with  a 

bar";  and 

9  -  Superscripts  in  parentheses  on  derivative  vectors?  elements  denote  the 

order  of  the  derivative. 

10  -  Formula  numbers  with  an  asterisk  are  rewritten  in  Appendix  I  along  with 
tie  corresponding  formulas  of  reference  1  whose  numbers  are  given  in 
brackets,  ° 


The  Matrix-Vector  Formulation  of  the  Classical 
Constrained  Least-Squares  Curve  Fitting  Procedure 


We  will  be  fitting  an  N-th  order  polynomial  to  a  span  of  (M»-l)  data 
points  x  ,  0<m<M,  M>N,  equally  spaced  with  an  interval  b.  The  Mtfr  position 
of  the  span  will  contain  the  current  (real-time)  data  point.  For  clearer 
delineation  we  will  use  the  notation  of  a  truncated  Taylor  series  expanded 
about  the  point  m*0  instead  of  a  polynomial  with,  undetermined  coefficients. 
Accordingly,  we  write: 
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where : 


A 

x  -  Smoothed  estimates  of  the  *,  the  Mfl  values  of  input  data 
m  contained  within  the  span,  iKra<M; 
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-  Estimates  of  the  derivatives  (parameters),  0<n<N,  evaluated 
at  m*»0. 


A(n) 

Given  the  we  wish  to  determine  the  values  of  the  x£  '  which  will  minimize 
the  squared-error  sum  \fi  : 
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In  matrix-vector  notation  (1)  and  (2)  are  written  as: 


A  A  .. 


x-e  «  x  »  Axr 


and 


T  *  T  A 

ip  »  e  e  »  (x-x)  (x-x) 
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Rearranging  (10),  substituting  (3),  and  solving  for  the  parameter  estimates 
gives  * 


or 


T  T  A 

A  x  -  ATAXq, 


^  -  [ATA]"1ATx, 


(ID* 


Substituting  (11)  into  (3)  we  obtain  the  desired  estimates  x  with  the  matrix- 
vector  formula  for  the  classical  unweighted  least-squares  curve  fitting 
procedure,  i.e., 


x 


AIaV^x. 


(12)* 


Now  suppose  we  have  the  very  practical  situation  that  we  require  not 
only  some  estimate  m  j  0.  within  the  span,  e.g. ,  the  mid-point,  but  also 
the  derivative  estimates  f  l<n<N.  is  given  by  the  Taylor  series  (1) 
after  the  parameter  estimates  xq  are  obtained  with  (11) .  The  derivative 
estimates  are  obtained  directly  by  successive  differentiation  of  (1),  i.e., 
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In  matrix-vector  notation  the  system  (13)  becomes 


/\ 


(14)* 
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d.  It  is  an  integration  operator. 


4>xr 


Substituting  (11)  into  (14)  gives  the  formula  for  evaluating  the  required 
estimates. 


m 


$m[ATA]"1ATx. 


(17)* 


The  foregoing  development  has  been  presented  to  demonstrate  our  matrix- 
vector  notation  in  obtaining  the  familiar  classical  least-squares  curve 
fitting  procedure.  The  constrained  procedure  begins  with  the  mathematical 
specification  of  the  "intercept"  and  "slope"  constraints.  The  intercept 
parameter  is  and.^the  slope  parameter  is  *q.  When  the  constraints  are 
not  applied, *q  and  xQ  are  evaluated  along  with  the  regaining  parameters 
using  (11) .  When  the  constraints  are  applied^  &q  and  *«  are  obtained  using 
(14)  with  mpl,  i.e., 


The  xq  are  the  parameter  estimates  obtained  in  the  preceding  computation  step. 
The  xq  are  the  predicted  values  of  the  current  parameter  estimates  which  are 
precisely  the  x^  obtained  in  the  preceding  computation  step.  We  apply  the 
intercept  and  slope  constraints  by  imposing  the  conditions: 
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toneidejationf  “17  ?artlCUl!,f  coa*“*1<>li*  is  justified  by  physical 
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specified  in  the  transition  matrix  $  in  (18).  F  C  1 

We  will  now  evaluate  the  remaining  parameter  estimates  unrfer  *ua 
constraint  conditions  «>.  Ihe  fit.c*,tep  is  to  p««iti“  7S  ^i„  (3) 
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It  Is  important  to  note  here  that  this  constrained  least-squares  curve¬ 
fitting  procedure  is  inherently  recursive  because  of  the  very  nature  of  the 
Intercept  and  slope  constraints,  i. a., using  information  obtained  in  the 
preceding  step  to  evaluate  the  estimates  in  the  current  step.  It  follows, 
then,  that  the  procedure  must  be  "initialized"  with  initial  constraint  values 
which  are  consistent  with  subsequent  data.  Any  appreciable  inconsistency 
generates  an  undesirable  transient  which  will  mask  out  the  desired  estimates 
until  the  transient  "settles  out".  In  practice,  the  initial  constraints  are 
obtained  with  an  unconstrained  curve  fit  over  the  first  M+l  data  points. 
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Matrix  -  Vector  Formulation  of  the  QD  Filter 


The  key  to  this  development  of  the  QD  filter  is  an  approximation 
requiring  reformulation  of  the  conditions  (28)  which  minimize  the  sum  of 
the  squared  errors  under  the  slope  and  intercept  constraints.  From  (28) 
we  have 


A2*c 


K2  W: 


(31)* 


From  (25)  and  (26)  we  can  write 


xc  ■  x-A1(x0)1-A2(x0)2+A2(^)2  ■  <Sx+A2(x0)2, 


where: 


6x  =  x-A1(xq)1-A2(x0)2. 


(33)* 


Substituting  (32)  into  (81)  and  rearranging  gives 


A^fix  -  A^AjSxq, 


(34)* 


where: 


6x0  5  <X0)2’(X0)2’ 


(35)* 


(x0)2  “  (xo>2+6xO* 


«WV*  >WW(  in-  «. 


Solving  (34)  for  the  fix^  gives 

6xq  -  tA2VV**. 


(37)* 


Substituting  (37)  into  (36)  gives 


(xQ)2 


(x0)2 


tA2TA2]“1A2T«x, 


(3ft)* 


in  which  the  unconstrained  parameter  estimates  (Xq)2  are  computed  as 
sums  of  the  predicted  values  and  corrections.  Clearly,  (38)  has  no 
computational  advantage  over  (29)  but  was  developed  to  provide  a  basis 
for  the  QD  approximation. 

We  observe  in  (33)  that  the  fix  vector  contains  the  last  (M*l)  input 
data  points,  whereas  the  QD  filter  uses  only  the  current  data  point.  If 
(38)  were  used  to  evaluate  the  unconstrained  parameter  estimates.,  an  x 
array  would  have  to  be  retained  to  evaluate  the  fix  in  each  computing  step. 

An  alternative  procedure  would  be  to  retain  a  fix  array  which  would  have  to 
be  shifted  in  each  computing  step  as  well  as  corrected  to  take  into  account 
the  new  values  of  parameter  estimates.  With  this  arrangement  each  input  data 
point  makes  a  contribution  through  fix  to  the  evaluation  of  the  parameter  estimates. 
Moreover,  each  data  point  makes  a  contribution  for  (M+l)  consecutive  computing 
steps  as  it  "passes  through  the  filter".  The  QD  alternative  is  to  let  each 
data  point  make  all  of  its  (Wfl)  contributions  during  the  current  computing 
step  and  drop  the  contributions  of  the  other  M  data  points.  To  implement  this 
reasoning  we  set  each  fix  element  equal  to  the  current  element  fix^  minus  the  sum 
of  the  corrections  it  would  receive  as  the  filter  processes  the  current  and 
subsequent  If  data  points.  For  example,  the  mth  element  of  fix  for  a  second- 
order  constrained  filter  would  be 


fix  =  fix  - 
m  M 
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(39) 
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in  which  the  fi>?0  .  is  the  correction  for  the  unconstrained  parameter  estimate 
used  to  compute  {he  correction  for  8xm  at  the  kth  computing  step  after  the 
corresponding  input  data  point  entered  the  filter. Since  each  ^  in  (39) 
can  not  be  determined  until  the  next  M  data  points- are  processed^  we  will  assume 
their  average  value  is  fix0t  the  parameter  correction  vector  which  is  being  evaluated 
during  the  current  computing  step.  Applying  this  assumption,  (39)  can  be  written 
in  the  matrix-vector  form; 
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In  order  to  obtain  an  explicit  expression  for  (6xft)«  in  this  modified  system, 

(43)  is  substituted  in  (34)  to  get  “ 

-  C(«^)2]  -  A2TA2  (6^)2  (46) 

A 

Solving  (46)  for  ($Xq)2  and  substituting  the  result  into  (36)  gives  the  QD 
unconstrained  parameter  estimates  (xq)2: 

(xq}2  "  (xq)2  +  IA2T(C+A2) ]"1A2Tv«xm.  (47) 

Using  (47)  in  (14)  gives  the  derivatives*  estimates  for  any  point  within  the  span: 


xm  “  *%  *  **  <*6>1 


(Xq)2 


(xq)2+[A2T(C+A2)  f1  A2Tv«xm 


[A2T(C+A^]"1A2Tv5xm 


where : 


A 

“  real-time  derivative  estimates  obtained  in  preceding 
computing  step. 


(53)* 


In  the  QD  filter  the  intercept  and  slope  constraints  are  applied  by  (52) 
in  which  <J>$  effects  only  the  corrections  on  the  unconstrained  parameters 
(xq)j  to  the  real-time  derivative  estimates.  (51)  is^written  in  several 
alternate  forms  to  show  that  the  parameter  estimates  Tq  evaluated  at  m*0 
in  the  preceding  step  are  implicitly  contained  in  the  OD  formulation. 


The  smoothed  QD  outputs  (i.e.,  those  nearest  the  point  of  constraint) 
are  given  by  (14) : 


A 


(54)* 


so  that  the  smoothed  QD  outputs  are  precisely  the  parameter  estimates  at 
the  point  of  constraint. 


Extensions  and  Applications  of  the  0D  Filter. 

It  is  a  very  simple  matter  to  mechanize  the  QD  filter  with  a  "variable 
span"  which  is  dynamically  self-ad.iusttng  (i.e.  .adaptive)  to  specified 
properties  of  the  current  data.  At  WSMR,  for  example,  the  QD  filter  in  the 
flight  safety  program  for  PERSHING  overflight  support  mentioned  earlier  uses 
about  a  three-to-one  variation  in  span  length  which  automatically  adjusts  to 
the  current  magnitude  of  the  total  acceleration,  i.e.»the  greater  the  accelera¬ 
tion,  the  shorter  the  span.  This  feature  effects  a  drastic  reduction  in  velocity 
overshoot  because  the  span  is  at  its  minimum  value  during  the  "staging"  of  this 
event.  Although  a  little  extra  storage  is  required  for  the  precomputed  coefficients 
for  each  span  length,  the  program  logic  for  this  feature  is  still  quite  simple. 
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It  should  be  pointed  out  that  variable  span  length  can  be  adapted  to  almost 
any  type  of  filter.  The  adaptation  to  QD  happens  to  be  extremely  simple. 

Reference  1  Includes  a  description  of  the  DFX  digital  filter  developed 
by  the  writer  in  the  Summer  of  1964.  Except  for  the  manner  In  which  the 
corrections  for  the  unconstrained  parameters  are  computed,  QD  and  DFX 
employ  the  same  basic  recursive  philosophy.  In  the  Spring  of  1965  the 
writer  presented  an  Informal  paper  to  the  Data  Reduction  and  Computer 
Working  Group  of  the  Inter-Range  Instrumentation  Group  at  the  Eastern  Test 
Range  (AFETR).  In  addition  to  explaining  the  basic  DFX  theory,  the  paper 
described  an  extremely  effective  technique  for  inserting  independently 
measured  derivative  information  into  the  filter  performing  its  normal 
filtering  functions.  More  specifically,  it  was  pointed  out  that  measured  \ 
velocity  data  could  be  simultaneously  inserted  into  a  DFX  filter  accepting' 
radar  posltlnn  data  to  essentially  eliminate  such  errors  as  "overshoot"  at 
motor  burnout  or  other  events  generating  abrupt  changes  in  acceleration.  \ 
To  do  this,  the  QD  or  DFX  predicted  velocity  estimate  ^  ie  replaced  directly 
with  the  corresponding  measured  velocity  component  v  so  that  the  predicted 
posltlnn  estimate  in  a  second  order  filter. 


*M  “  *M  +  T(\  +  \)» 


(55) 


would  be  replaced  by 


M 


(56) 


It  may  be  argued  that  this  procedure  does  not  buy  anything  because  the  velocity 
data  are  already  available.  This  is  an  erroneous  conclusion  because,  with 
Inserted  velocity  data,  the  filter  span  cah  be  significantly  Increased  to 
provide  dramatic  Improvement  in  filtering  the  radar  data  with  no  Increased 
distortion  of  the  true  information  in  the  radar  data.  In  effect,  then,  the 
filter  is  operating  on  the  difference  between  the  trajectory  described  by 
the  radar  data  and  the  trajectory  described  by  th^  integration  of  the  velocity 
data.  This  unique  capability  has  not  been  used  in  practice  at  WSMR  because 
suitable  measured  velocity  data  have  not  .been  available.  The  writer,  however, 
has  tested  this  technique  in  a  wide  variety  of  simulated  situations  with 
convincing  results.  This  same  improvement  in  performance  has  also  been  noted 
when  velocity  data  along  with  position  data  have  been  Inserted  into  the 
Kalman  filter. 
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At  the  AFETR  meeting  in  the  Spring  of  1965,  the  writer  also  described 
how  this  technique  could  be  extended  to  the  insertion  of  acceleration  data. 

In  filtering  missile  trajectory  data,  for  example,  measured  acceleration 
data  “X  obtained  via  telemetry  can  be  Inserted  into  the  filter  along  with 
the  measured  position  data  by  replacing  the  predicted  velocity  formula 

\  •  \  +  «7) 


with 


\  -  Vhv  <58> 


A  special  case  of  this  acceleration  data  insertion  has  been  used  very 
successfully  at  WSMR  for  more  than  two  years  in  the  flight-safety  support 
program  mentioned  earlier.  In  this  case,  a  burnout  signal  generated  from 
telemetered- data  is  made  available  to  the  real-time  computer  program. 

When  this  burnout  signal  arrives  the  Tty  components  of  the  QD  filters  are 
replaced  by  computed  gravity  components  for  the  next  several  computing 
cycles.  In  other  words,  when  the  burnout  event  occurs,  the  QD  filters 
are  forced  to  forget  about  acceleration  history.  This  procedure  essentially 
eliminates  the  Instantaneous  impact  prediction  overshoot  problem.  This 
procedure  also  suggests  that  much  more  effective  filtering  of  radar  data  can  be 
obtained  from  a  ballistic  missile  with  negligible  drag  by  inserting  the 
gravity  components  into  the  QD  filters. 

It  should  be  pointed  out  here  that,  when  velocity  or  acceleration 
data  are  inserted  into  a  QD  filter,  the  smooth  data  formula  (54)  is  no 
longer  valid  because  the  -Inserted  data  are  arbitrary  as  for  as  the  filter 
is  concerned.  (54)  is  based  on  a  constant  N-th  order  derivative  across 
the  filter  span, which  condition  is  obviously  violated  with  the  insertion  of 
arbitrary  derivative  data. 

We  have  also  experimented  extensively  with  QD  filter  arrangements  using 
a  simple  deterministic  mathematical  model  to  generate  the  Inserted  acceleration 
data.  As  in  the  Kalman  filter,  the  effectiveness  of  this  arrangement  is  highly 
dependent  on  how  well  the  trajectory  generated  by  the  model  matches  the  measur¬ 
ed  trajectory.  Feedback  techniques  have  been  used  to  force  the  measured 
trajectory  data  to  dynamically  update  (correct)  the  model  coefficient*.  Includ¬ 
ing  a  model  with  the  QD  filter  can  provide  performance  comparable  to  the  Kalman 
filter  in  relatively  simple  applications,  e.g.,«m  N-station  radar  data  reduction. 
When  the  application  becomes  complek,  e.g.,the  dozens  of  coefficients  contained 


18 


cowputing  efficiency  ““the  QD%u£^th‘iMMMd“UlMtlve  d^'^l]'''  th' 
and  is  certainly  beywd  the  scopt  of  this'report."0"''7  °'  eitten8lv'  Investigation 


A 

X  ■ 

m 

.nT 

$  Xq 

(14) 

[10J 

:<j) . 
xi 

a  f  (ih)n-J  "(n) 

n-j  <^>l  *»  *  J“°.  1.  •  -N. 

* 

*  *  *  * 

A 

X  = 

m 

$m[ATA]"'1ATx 

(17) 

[11] 

:<j) 

Xi 

?  (1) 

k=0  b;Lk  Sk’  ^  "  °»  X»  *  *  *  N»  1  ■  !»  2»  -  .  .  M 

[12] 

b«> 

Dik 

„  y  (ih)n"j  r 
n-j  ^n_ j  )  1  nk 

* 

*  *  *  * 

d\jj 

W2 

A  _  A 

*  -2(xc-xc)  A2  *  “2A2(VXc)  “  0 

(28) 

[15] 

?  1  A 

l  ,m  <VV  i  -  2,  3,  ...  N. 

in— x 

* 

*  *  *  * 

A 

(ltO>2 

*  [A^f1  A^c 

(29) 

[19] 

;<j>. 

xo 

N 

L  3  - 2- 3 —  n-- 

[18] 

V 

I  J  l  mj  -hxQ  l  mi+l 

ra3!  m»l  jn-i 

22 


6x0  “  <xO>2"(xO)2 


6jf0j> 

-  4  „7  , 
*0  Xq  •  J  2»  3. 

* 

*  * 

6x0  - 

sSJ) 

N 

■  La^ 

V 

M 

I  m3  Ax  ,  j.  2,  3,  . 
m-1 

* 

*  * 

A 

(V2 

-  (^j)2  +  J'1AT26x 

C<j) , 

*m>+  1$ Wv  J- 

criK 

^2  (n-j)!  V  *  "  0 

r(J)  . 
Sk 

y  (roon'3  - 

£«j  (n-j)!  V  J  “  2 

* 

*  * 

6x  =  v6xw  - 


M  ~  C(5V^ 


Axi  ■  %-  l  I  9?  ««n!J 


A 


*M  “ 

xM  +  w4>k 

(50) 

[68] 

A(j) 

V  1 

■  ^3>  +  «j  *V  J  -  0,  1,  .  . 

.  N 

[69] 

*1  - 

L  (^j)t  V  P  «  2  for  j  - 
n-p  J  P  -  j  for  j  - 

o,  1 

2,  3,  ...  N 

* 

[66] 

QJ  - 

N  H 

L  qj"  L  3  ■  2*  3>  • 

.  .  N 

<v 

■  'V1 

[64] 

qJn  ’ 

hn  M  M 

■  5r  I,”3  l  p" 

o-l  p-m 

* 

*  *  * 

* 

*Ha 

4  - 

(51) 

[30] 

N  .  n-j  .  » 

n-j  **  "  °*  X»  *  * 

.  N 

* 

*  *  * 

* 

w  » 

r 

$2  [A^CCfAj)]"1  a£v 

(52) 

[69] 

^  • 

f  (Mh)n"jn  *  —  _ 

*  (K^j)!  V  p  "  2  for  J  *  °» 

n*P  J  P  -  j  for  j  -  2, 

1 

3t  •  •  •  N 

[66], 

[64]  above 

* 

*  *  * 

* 

25 


J‘JM!l.«  A&  AM  J  liifiltfi  W  *  JiflJifiuSsK 


References 


T“m,r 


i 

» 

l 


1.  "QD  -  A  New  Efficient  Digital  Filter”,  W.  A.  McCool,  Internal  Memorandum 
No.  60,  Aug.  1967,  Analysis  and  Computation  Directorate,  White  Sands  Missile 
Range,  New  Mexico. 

2.  "Optimizing  the  Dynamic  Parameters  of  a  Track-While-Scan  System”,.  Jack 
Sklansky,  RCA  Review,  June  1957,  RCA  Laboratories,  Princeton,  N.J. 

3.  "Synthesis  of  an  Optimal  Set  of  Radar  Track-While-Scan  Smoothing  Equations", 
T.R.  Bendict  and  G.  W.  Bordner,  IRE  Transactions  on  Automatic  Control,  1962. 

4.  "Development  of  a  Linear  Recursive  Filter" *  J.  J.  Lynn,  Oct.  1964,  In-Flight 
Analysis,  Mathematical  Services,  RCA  Service  Co. 


fi 


3 


.1 

-4 


27 


